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ABSTRACT 

Given a set of images, whose pixel values can be considered as the components of a 
vector, it is interesting to estimate the modulus of such a vector in some localised 
areas corresponding to a compact signal. For instance, the detection/estimation of 
a polarized signal in compact sources immersed in a background is relevant in some 
fields like astrophysics. We develop two different techniques, one based on the Neyman- 
Pearson lemma, the Neyman-Pearson filter (NPF), and another based on prefiltering- 
before-fusion, the filtered fusion (FF), to deal with the problem of detection of the 
source and estimation of the polarization given two or three images corresponding to 
the different components of polarization (two for linear polarization, three including 
circular polarization). For the case of linear polarization, we have performed numerical 
simulations on two-dimensional patches to test these filters following two different 
approaches (a blind and a non-blind detection), considering extragalactic point sources 
immersed in cosmic microwave background (CMB) and non-stationary noise with the 
conditions of the 70 GHz Planck channel. The FF outperforms the NPF, especially 
for low fluxes. We can detect with the FF extragalactic sources in a high noise zone 
with fluxes ^ (0.42,0.36) Jy for (blind/non-blind) detection and in a low noise zone 
with fluxes ^ (0.22,0.18) Jy for (blind/non-blind) detection with low errors in the 
estimated flux and position. 

Key words: polarization - methods: data analysis - techniques: image processing - 
cosmic microwave background - radio continuum: galaxies. 



1 INTRODUCTION 

The detection and estimation of the intensity of compact 
objects (or small regions) embedded in a background plus 
instrumental noise is relevant in different contexts, e.g. as- 
trophysics, cosmology, medicine, teledetection, radar, etc. 
Different techniques have proven useful in the literature. 
Some of the proposed te chniques are f r equentist, such as the 
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L2009) that correspond to scalar, vector or matrix filters, re- 
spectively. Other frequentist techniques include continuous 



wavelets like the standard Mexican Hat ijVielva et al. I I2OOII , 
l2003l : lBarnard et al . "-'I HI |':'s anz e t al]|2006l) an d other mem- 
bers of its family (jGonzalez-Nu cvo et al.l [20061 and, more 
generally, filters based on the Neym an-Pearson approach 
using the distribution of maxima l|L6pez-Caniego et al.l 
l2005al : iLopez-Caniego et all l2005bD . All of these filters 
have been used in the literature, in particular, for 
the detection and estimation of the intensity of point- 
like sources (i.e. extragalactic objects) in Cosmic Mi- 
crowave Background (CMB) maps. In addition, some of 
them have been applied to real data like those obt ained 
by the WMAP satellite (Lopez-Caniego et al.' '200 7|) and 
simu lated data (Lopez-Canieg o et al.. ,2006 : Leach et al.l 
l2008h for the upcomin g experiment on board the Planck 
satellite (|Tauberl 120051 ). Besi des, Bayesian methods have 
also been recently d eveloped (|Hobson fc McLachlanI l2003l : 
ICarvalho et~aLll2009l ). 
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In some applications it is important to measure not 
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only the intensity of the light (signal) but also its polar- 
ization. An example is the study of CMB radiation. The 
polarization is given by the Stokes parameters Q,U,V , and 
the tot al intensity of polarizatio n is P = (Q^ + [/^ + 
(|Kamionkowski et al.|[l997h . For linear polarization, 
the previous expression reduces to P = {Q'^+U'^y^^ . In such 
cases, three or two images are added quadratically followed 
by a square root. 

Let us consider the case of linear polarization. In this 
case we have two images Q, U and different approaches 
can be used to deal with detection/estimation of point-like 
sources on these maps. On the one hand, one can try to get 
the polarization P directly on the P-map. In this approach, 
we will consider one filter, obtained through the Neyman- 
Pearson technique (NPF). On the other hand, we can oper- 
ate with two matched filters, each one on Q and U followed 
by a quadratic fusion and square root. We will call this pro- 
cedure filtered fusion (FF). It is clear that from a formal 
point of view we are trying to ask which is the optimal way 
to find the modulus of a vector given the components. In 
the case we have only the map of the modulus of a vector 
and the components are unknown, the FF cannot be applied 
and the only possibility is the NPF. 

We will develop the methodology for the cases of a 2- 
vector and a 3- vector because of the possible interesting ap- 
plications to the 2-plane and 3-space . We will show the re- 
sults when using numerical simulations on flat patches that 
are relevant for the component separation of sources (linear- 
polarized extragalactic sources) in CMB maps. 

In section [21 we develop the methodology. In section |31 
we describe the numerical simulations done in order to test 
the previous techniques. In section 2] we present the main 
results and in section [S] we give the main conclusions. 



2.1.1 Neyman-Pearson filter (NPF) on the P-map 

If the noise is distributed normally and independently, then 
at any point the integration over the polar angle leads 
to the 2D-Rayleigh distribution of P in absence of the 
source (|Papouli s 1984) 



/(PIO) = -.e 



(2) 



whereas if the source is pre sent , with amplitude A, one ob- 
tains the Rice distribution (|Ricelll954l ) 



j{P\A) = ^e- 



'44) 



(3) 



where Jo is the modified Bessel function of zero order. If our 
image is pixelised, the different data Pi, i = 1, . . . , n, with 
n the number of pixels, will follow the two distributions 
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being cr; the value of a in the i' pixel, i^o, Hi the null 
hypothesis (absence of source) and the alternative (presence 
of source), respectively, and n the profile at the i"' pixel. 
The log-likelihood is defined by 



UA|P) = log 
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The maximum likelihood estimator of the amplitude. A, is 
given by the solution of the equation 

= —jt- (7) 
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2 METHODOLOGY 



This equation can be interpreted as a non-linear filter op- 
erating on the data that we shall call the Neyman-Pearson 
filter (NPF). 



2.1 2- vector 

To develop our methodology, we shall assume that we have 
a compact source, located for simplicity at the centre of two 
images Q, U and characterised by amplitudes Aq, Au and 
a profile r{x), immersed in noises nQ{x) and nu{x) that 
are Gaussian and independently distributed with zero mean 
and dispersions o-q{x) and au{x). Again, for simplicity, we 
will consider that o-q{x) — au{x) — cr{x), a condition that 
is verified in most polarization detectors. In general, we will 
consider that the noise is non-stationary. We remark that the 
previous assumptions can be easily generalised to different 
profiles and different types of noise in the two images but 
we will not consider it in this paper. We will assume a linear 
model for the two images 



dq^uix) = Aq,ut(x) + nQ,u{x). 



(1) 



The P-map, P(f) = {Q^[x) + U'^[x)f^'^, is charac- 
terised by a source at the centre of the image with amplitude 
A = {Aq + A^Y^^ immersed in non-additive noise which is 
correlated with the signal. 



2.1.2 Filtered fusion (FF) 

In this case we use the same matched filter (MF) oper- 
ati ng over the two imag es Q, U, respectively, as given 
by lArgiieso fc Sanj|200S 
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Then, with the two filtered images Qmf, Umf we make the 
non- linear fusion P = {Q%ip + Ump)^^^ pixel by pixel. 

2.2 3- vector 

Now, we shall assume that we have a compact source at the 
centre of three images Q, U, V characterised by amplitudes 
Aq, Au, Av and a profile t{x) immersed in noise nq^ij^vix) 
that is Gaussian and independently distributed with zero 
mean and dispersion a{x). In general, we will consider that 
the noise is non-stationary. We will assume a linear model 
for the three images 



• ,u,v{x) = AQ^u,vr{x) + nQ^u,v{x). 



(9) 
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The P-map, P(x) = (Q^ix) + U^{x) + y2(f))^''^ is charac- 
terised by a source at the centre of the image with ampUtude 
A = {A'^ + A^ + AY)^^^ immersed in non-additive noise cor- 
related with the signal. 



2.2.1 Neyman-Pearson filter (NPF) on the P-map 

If the noise is distributed normally and independently, then 
at any point the integration over the angles leads to the 3D- 
Rayleigh distribution in absence of the source, also called 
the Maxwell-Boltzmann distribution in Physics 
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whereas if the source is present, with amplitude A, one ob- 
tains the distribution 



If our image is pixelised, the different data Pi, i = 1, . . . ,n 
follow the two distributions 
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being Hq, Hi the null hypothesis (absence of source) and 
the alternative (presence of source), respectively and n the 
profile at the i**^ pixel. The log-likelihood is defined by 



-nlogyl + log 
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(14) 



The maximum likelihood estimator of the amplitude. A, is 
given by the solution of the equation 



^yiCot\i{Ayi), yi = ^-^. (15) 



This equation can be interpreted as a non-linear filter oper- 
ating on the data that we shall also call the Neyman-Pearson 
fiher (NPF). 

2.2.2 Filtered fusion (FF) 

In this case we use the same MF operating over the three 
images Q, U, V, respectively, as given by equation ([8]). Then, 
with the three filtered images Qmf, Umf, Vmf we make the 
non-linear fusion P = {Q%iF + U^iF + '^MFy^^ pixel by pixel. 



3 SIMULATIONS 



The E uropean Space Agency (ESA) Planck satellite (jTaubed 
I2OO5I ). to be launched in 2009, will measure the anisotropies 
of the CMB with unprecedented accuracy and angular reso- 
lution. It will also analyse the polarization of this primordial 
light. It is of great interest the detection and estimation of 
polarized sources in CMB maps (|The Planck Collaboration! 



l2006l : lTucci et al.ll2004 120051): since t his radiation is linearly 
polarized I Kamionkowski et al.lllQOTl '), V = 0, we will apply 
the methods for the detection/estimation of the modulus of 
a 2- vector presented in section [2TT] However, some cosmolog- 
ical models predict a pos sible circular polarization of CMB 
radiation (see for example lCoorav et al.l[20o3 : lAgarwal et al.l 
I2OO8I I. Even if CMB is not circularly polarized, the extra- 
galac tic radio sources can indeed show cir c ular polariza- 
tion (jAller et al.l l2005l : iKirk fc Tsa"n3 I2OO6I : iHoman et aD 
l2006h . Besides, circular polarization occur s in many other 
astrophys ical areas, fro m Solar Physics (iTritschler et al 



2007: RcinereEaD[2003) to interestellar medium (ICox et al 



2003), just to put a few examples. In all these cases the re- 



sults for the modulus of a 3- vector presented in section [2^ 
could be useful. 

In order to compare and evaluate the performance of 
the filters presented in section [2m we have carried out two- 
dimensional simulati ons with the characteristics of t he 70 
GHz Planck channel l|The Planck Collaborationll2006l ). The 
simulated images have 16 x 16 pixels with a pixel size of 3.43 
arcmin. We simulate the Q and U components of the linear 
polarization as follows: each component consists of Gaus- 
sian uncorrelated detector noise plus the contribution of the 
CMB and a polarized point source filtered with a Gaussian- 
shaped beam whose full width at half maximum (FWHM) 
is 14', (the FWHM of the 70 GHz Planck channel beam). 
So the source polarization components can be written as 



SQ = Aq exp 
su = Au exp 



m 

272 
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272 



(16) 
(17) 



where 7 is the beam dispersion (size) and we assume that the 
source is centred at the origin. The CMB simulation is based 
on the observed WMAP low multipoles, and on a Gaussian 
realisation assuming the WMAP best-fit Ce at high multi- 
poles. We do not include other possible foregrounds, since we 
are doing a first approximation to the detection/estimation 
problem and we assume that we apply our filters to rela- 
tively clean areas far away enough from the Galactic plane. 
Alternatively, we could consider a case in which the major- 
ity of foregrounds have been previously removed by means 
of some component separation method. 

We consider realistic non-stationary detector noise. We 
have simulated the noise sky pattern for a Planck flight du- 
ration of 14 months, assuming a simple cycloidal scanning 
strategy with a 7 degree slow variation in the Ecliptic coUat- 
itude of the z-axis. This scanning strategy implies that the 
sky will be covered inhomogeneously. The si mulations have 
t he sa me characteristics as the ones used in lArgiieso et al.l 
()2006l ). In order to illustrate the effects of non-stationary 
noise we have chosen two representative zones of the sky: 
one zone of high noise but quite isotropic and another zone 
of low noise but more anisotropic. The first zone is located 
in a region far from the Ecliptic poles, where the noise pat- 
tern is quite uniform and the number of hits per pixel of 
the detector is small. The average r.m.s deviation of the 
first zone (high noise) in units of AT/T (thermodynamic) 
is (J = 3 X 10~^ (for each component Q and U) and its 
standard deviation is 4.2 x 10~®. For this particular scan- 
ning strategy approximately 25% of the sky has this kind of 
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noise pattern. The second zone is close to one of the Ecliptic 
poles, where the sky is scanned more times (low noise level) 
but the hit pattern is very inhomogeneous. It has an average 
r.m.s. deviation a — 1.1 x 10~^ and its standard deviation is 
3.8 X 10"^. Then, in the second zone the noise is lower but 
proportionally more anisotropic. For this scanning strategy, 
~ 6% of the sky has a noise pattern with these character- 
istics. Other zones of the sky would be intermediate cases 
between those considered here. 

We take values of the source fluxes (before filtering with 
the Gaussian beam) Fq and Fu ranging from 0.1 Jy to 0.5 
Jy with a step of 0.1 Jy. A flux of 0.1 Jy corresponds for 
the 70 GHz channel to AT/T = 1.9 x 10"^ so that it is of 
the order of the detector noise r.m.s deviation. The fluxes 
also corres pond to typical polarization fluxes (|Tucci et al.l 
I2OO4I . [2OO5I ). As we will see in section |4j for some low-noise 
cases it is necessary to simulate even lower flux sources in 
order to study the behaviour of the filters in the low signal 
to noise ratio regime. We will explicitly report on this in the 
appropriate cases in section |4] The number of simulations is 
500 for each combination of pairs of values of Fq and Fu. 
After carrying out the corresponding simulations for Q and 
U, we add them quadrati cally and take the square root to 
calculate P = ^ + , the polarization modulus. 

We study two different detection types: blind detection 
and non-blind detection. In the former case, we assume that 
we do not know the position of the source and then we place 
it at random in the image, in the latter case we know the 
source position and then we place it at the centre of the 
patch. We have considered images of 16 x 16 pixels in order 
to do fast calculations. In order to avoid border effects, we 
simulate and filter 24 x 24 pixel patches and, after the fil- 
tering step, we retain only the 16 x 16 pixel central square. 
We have tested the case of larger patches but the obtained 
results do not change significantly. 

We use two different filters: the filtered fusion (FF), 
which consists in the application of the matched filter to the 
images in Q and U separately in a first moment, and then the 
calculation/construction of P from the matched-filtered im- 
ages Qmf and Umf and the Neyman- Pearson filter (NPF), 
applied directly on P, derived from the Neyman-Pearson 
lemma and presented with detail in section [2] These filters 
are suitable for the case of uncorrelated noise, but we apply 
them to simulations including the CMB, which is correlated. 
However, we have checked that our results are similar when 
we consider simulations with and without the CMB: the rel- 
ative differences of the errors in the estimated fluxes and 
positions are at most of a few percent. This is due to the 
low value of the CMB r.m.s. deviation cr ~ 6 x 10~^ as 
compared to that of the detector noise. Hence, the meth- 
ods derived in section 12.11 are also suitable for simulations 
including the CMB in the 70 GHz Planck case we consider. 

In the blind case we apply these filters to each simula- 
tion, centering the filters successively at each pixel, since we 
do not know the source position. We estimate the source am- 
plitude A for the NPF, in this case we calculate the value 
of A which maximises the log- likelihood, eqs. ^ and ([7]). 
For the FF we estimate separately Aq and Au and obtain 

a = ,/a^Ta^. 

We compute the absolute maximum of A in each filtered 
map and keep this value as the estimated value of the po- 
larization P of the source and the position of the maximum 



as the position of the source. Note that in the more realis- 
tic case where more than one source can be present in the 
images, it is still possible to proceed as described by look- 
ing for local peaks in the image. In the non-blind case we 
only centre the filters in the pixels included in one FWHM 
of the source centre (approximately 10% of the total). In 
this way, we use the knowledge of the source position, then 
we calculate the absolute maximum of the estimated A in 
these pixels. We also calculate the significance level of each 
detection. In order to do this, we carry out 1000 simulations 
with Aq = and Au = 0, and we calculate the estimated 
value of the source polarization in this case for each filter. 
We consider the null hypothesis, Hq, there is no polarized 
source, against the alternative hypothesis, Hi, there is a 
polarized source. We set a significance level a — 0.05, this 
means that we reject the null hypothesis when a simulation 
has an estimated source amplitude higher than that of 95% 
of the simulations without polarized source. We define the 
power of the test as 1 — 5, with 5 the probability of accept- 
ing the null hypothesis when it is false, i.e. the power is the 
proportion of simulations with polarized source with an esti- 
mated amplitude higher than that of 95% of the simulations 
without source. The higher the power the more efficient the 
filter is for detection. Note that the test can be performed 
in the same way in the blind and non-blind cases: in the 
second case we know the position of the source, but we do 
not know whether it is polarized or not. 

4 RESULTS 
4.1 Blind case 

We carry out simulations in the blind case for the high noise 
and low noise zones as explained above. We apply the filters 
to the images and calculate the absolute maximum of A 
for each filtered image. If the detection has a significance 
higher than a = 0.05 we consider it as a real detection, 
otherwise it would be a spurious one. For this significance 
level, we calculate the power of the detection test for the 
different pairs of Fq and Fu values. We also calculate the 
estimated value of the polarization amplitude A, convert it 
to the estimated ffux in Jy, F, and compare it with the 
real value F = \/P'q + ^u- P^*^^ ^ against F for the 
NPF and the FF with high and low noise in Fig. [T] We also 
compute the relative error of this estimation and its absolute 
value. For each pair of Fq and Fu values, we obtain the 
average and 68% confidence intervals of the absolute value 
of the relative error, taking into account the 500 performed 
simulations. These values are plotted against F in Fig. [2] 
for the same cases shown in Fig. [T] We also calculate the 
estimated position of the detected source and obtain the 
position error expressed in terms of the number of pixels. 
The power and the average of the relative error in F, of 
its absolute value and of the position error are presented 
in Table [1] for the two different types of noise. The rows in 
the Table are sorted in ascending order of F. For the high 
noise case, we see in the Table that the power for the FF 
is higher than for the NPF. The position and polarization 
errors are also lower for the FF, this can also be seen in 
Figs. [1] and [1 For F > 0.6 Jy, the FF and NPF perform 
similarly. For ffuxes higher than 0.42 Jy, the power with the 
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Figure 1. Blind detection. Estimated flux F in Jy of the polarized sources plotted against their real flux F. The average and 68% 
confidence intervals (vertical bar) of 500 simulations are plotted. Top left: the NPF has been used and the noise corresponds to the high 
noise zone. Top right: FF and high noise. Bottom left: NPF and low noise . Bottom right: FF and low noise. In all the plots the straight 
line F = F is drawn for comparison. 



FF is = 100%, the average relative error (bias) is ^ 0.03, the 
average of the absolute value of the relative error is 0.12 
and the average position error is 5C 0.36. From now on, we 
will use the flux limit for which the power is 100% as a 
measure of the filter performance. 

The results for the low noise case are logically much bet- 
ter, the FF also performs better than the NPF. The power 
of the two filters quickly reaches 100% for fluxes F > 0.28 
Jy and the errors in both flux and position remain stable 
from fluxes F 0.40 Jy. We therefore have cut the Table 
at = 0.40 Jy. In order to have a better sampling of the 
interesting signal to noise regime, we have simulated (using 
the same number of simulations as in the other cases) in 
addition the flux pairs {Fq,Fu) = (0.05,0.10), (0.00,0.15), 
(0.05, 0.15), (0.10, 0.15) and (0.00, 0.25) Jy. This way the Ta- 
ble is much more informative. For fluxes higher than 0.22 Jy, 
the power with the FF is 100%, the average bias is ^ 0.01, 
the average of the absolute value of the relative error is 
0.09 and the average position error is ^ 0.21. 

The FF is also much faster than the NPF. For instance 
the analysis of an image of 64x64 pixels in a personal com- 
puter takes 7 seconds with the FF and about 6 minutes 
with the NPF. This is due to the maximization involved. 
The computation time grows proportionally to the number 
of pixels, so that the NPF could be too slow if we want to 
analyze large images. This is the reason why we have con- 
sidered small patches. 



Finally, we show in Fig. |3] four images corresponding to 
a polarized source with {Fq,Fu) — (0.4,0.4) embedded in 
high noise. For the sake of a better visualisation, we show 
64 X 64 pixel images instead of the 16 x 16 sized images 
used in the simulations. We show the original image in P 
including noise, CMB and source, the image of the source, 
the image filtered with the NPF and the image treated with 
the FF method. By simple visual inspection, it is possible 
to see that the the FF performs better than the NPF. 



4.2 Non-blind case 

We carry out simulations with a polarized source placed in 
the image centre and filter the images with the two different 
filters, centering them in the pixels at a distance from the 
source less than one FWHM. We simulate the same range 
of [Fq,Fu) pairs as in the blind case (section 14. ip . Then, 
as in the blind case, we calculate the maximum in A, using 
only the selected pixels. This maximum is the estimated po- 
larization of the source A and the pixel where we find the 
maximum is the source position, we convert this amplitude 
to the estimated flux in Jy, F. We have also performed a 
detection test, accepting as real sources only those detected 
with a significance lower than a = 0.05. In Tabled we write 
the power of the test, the average of the relative polariza- 
tion error, of its absolute value and of the position error. Our 
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Figure 2. Blind detection. Absolute value of the relative error, |err|, in the estimation of the flux of polarized sources plotted against 
the real fiux F in Jy. The average and 68% confidence intervals (vertical bar) of 500 simulations are plotted. Top left: The NPF is used 
with high noise. Top right: FF and high noise. Bottom left: NPF and low noise. Bottom right: FF and low noise. 



results are obtained from 500 simulations for each combina- 
tion of pairs {Fq , Fu) with high and low noise. The detection 
power is higher for the FF than for the other filter. For the 
high noise case, the improvement is very clear for F ^ 0.50 
Jy; we also obtain higher powers for the non-blind case than 
for the blind one. This could be expected, since we know the 
source position in the former case. 

The position and polarization errors are also lower for 
the FF, this can also be seen in Figs. |4] and O which show 
the same quantities as Figs. [T] and [5] for non-blind detection. 
For fluxes higher than 0.36 Jy, the power with the FF is 
100%, the bias is ^ 0.03,the average of the absolute value 
of the relative error is ^ 0.13 and the average position error 
is ^ 0.55. The flux and position errors are lower in the non- 
blind case than in the blind one, especially for low fluxes. 

The results for the low noise case are much better, the 
FF outperforms the NPF. From 0.28 Jy on the results are 
quite similar. For fluxes higher than 0.18 Jy, the power with 
the FF is 100%, the average bias is 5C 0.02, the average 
of the absolute value of the relative error is 5C 0.11 and the 
average position error is ^5 0.39. Similarly to what happened 
with the low noise case of Table [1] we have cut the Table 
at flux F = 0.40 Jy and we have added new cases in order 
to have a better sampling of the interesting signal to noise 
regime. In this case, we have simulated the pairs {Fq,Fu) = 
(0.05,0.10), (0.00,0.12), (0.00,0.13), (0.00,0.15), (0.10,0.15) 
and (0.00, 0.25) Jy. The non-blind case gives better results 
than the blind one, especially for low fluxes. 



5 CONCLUSIONS 

In this paper, we deal with the detection and estimation 
of the modulus of a vector, a problem of great interest 
in general and in particular in astrophysics when we con- 
sider the polarization of the cosmic microwave background 
(CMB). The polarization intensity P is defined as P = 
{Q^ + U^ + V'^)2, where Q, U and V are the Stokes parame- 
ters. We consider the case of images in Q, U and V consisting 
of a compact source with a profile t{x) immersed in Gaus- 
sian uncorrelated noise. We intend to detect the source and 
estimate its amplitude in P by using two different methods: 
a Neyman-Pearson filter (NPF) operating in P and based 
on the maximisation of the corresponding log-likelihood and 
a filtered fusion (FF) procedure, i.e. the application of the 
MF on the images of Q, U and V and the combination of 
the corresponding estimates by making the non-linear fusion 
P = {Q\ip 4- Ulip + Vmp)2 . We present the two filters in 
section[2]for two-dimensional, = 0, and three-dimensional 
vectors, deriving the corresponding expressions for the esti- 
mation of the polarized source amplitude. 

Since we are interested in applying these different fil- 
ters to the CMB and this radiation is linearly polarized, we 
will only consider the two-dimensional vector case in our 
simulations. 

Our goal is to compare the performance of the filters 
when applied to simulated images. Then, we carry out two- 
dimensional simulations with the characteristics of the 70 
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Figure 3. Blind detection image. Top left: Image of a polarized source filtered with a Gaussian beam (FWHM=f 4') placed in the centre 
of an image of 64 X 64 pixels with pixel size=3.43'. The source polarization flux components are {Fq, Fjj) = (0.4,0.4) in Jy and it is 
embedded in CMB plus detector noise (high noise zone). Top right: The image of the polarized source only. Bottom left: The first image 
after the application of the NPF. Bottom right: The first image after the application of the FF. 



GHz Planck channel. The images have 16 x 16 pixels with a 
pixel size of 3.43'. We simulate the Q and U components con- 
sisting of a compact source filtered with a Gaussian-shaped 
beam (FWHM of 14') plus CMB and non-stationary detec- 
tor noise. We consider two typical zones of the Planck sur- 
vey: one with high noise and quite isotropic and another one 
with low noise but proportionally more anisotropic. These 
zones are extreme cases for the assumed scanning strategy 
we have chosen, and any other zone of the sky is an inter- 
mediate case between these two. 

We study two types of detection: blind detection, in 
which we do not know the source position and non-blind, 
in which the position is known; we place the source at the 
centre of the image. We take values of the source fluxes in 
Q and U, Fq, and Fu, ranging from 0.1 Jy to 0.5 Jy with 
a step of 0.1 Jy. Note that for extragalactic objects both Q 
and U can take negative values, but since the sign of both 
components is irrelevant for the calculation of P here we only 
simulate the positive case. We carry out 500 simulations for 
each pair {Fq,Fu) and for all the cases of blind and non- 
blind detection with high and low noise. 

We apply the fllters to the simulated images, estimat- 
ing the source amplitude A as the maximum value of P in 
the filtered images and the source position as the position of 
this maximum, then we convert the source amplitude to the 
source flux in Jy . We fix a significance a — 0.05 for the de- 
tection and calculate a detection power for this significance. 



see section [3] for more details. We also calculate the relative 
error of the estimated flux, its absolute value and the posi- 
tion error (in number of pixels), these errors together with 
the detection power are written in Tables [T] and [5] for the 
blind and non-blind case. We also show the estimated fluxes 
and the absolute value of the relative errors in Figs. [T] (2] ID 
and [S] for the different cases. 

The FF performs better than the NPF (as can be seen 
in the Tables and Figures), especially for low fluxes and 
it is also much faster than the NPF. However NPF is still 
interesting in a case where only the modulus of a vector is 
known and not its components (for example, if we had a map 
of P polarization but not the Q and U maps separately) . The 
powers are much higher and the errors much lower for the 
low noise case than for the high noise one. The filters also 
perform better in the non-blind case than in the blind one, 
especially for low fluxes. 

We can detect extragalactic point sources in polariza- 
tion images (at 100% power) with the FF in the high noise 
zone with fluxes ^ (0.42,0.36) Jy for (blind/non-blind) de- 
tection and in the low noise zone with fluxes ^ (0.22, 0.18) 
Jy for (blind/non-blind) detection. The bias and the posi- 
tion error are very low ^ 0.03 and < 1 pixel, respectively 
for all these fluxes. 
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Figure 4. Non-blind detection. Estimated flux of the polarized sources _F in Jy plotted against their real flux F. The average and 68% 
confidence intervals (vertical bar) of 500 simulations are plotted. Top left: the NPF has been used and the noise corresponds to the high 
noise zone. Top right: FF and high noise. Bottom left: NPF and low noise . Bottom right: FF and low noise. In all the plots the straight 
Une F = F is drawn for comparison. 
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Table 1. BLIND DETECTION. First column: pairs of values of Fq and Fu in Jy and the corresponding value of F for the high-noise 
zone (top) and the low-noise zone (bottom). Second and third columns: detection power for the NPF and FF (percentage). Fourth 
and fifth columns: flux relative errors (average from 500 simulations) for the two filters. Sixth and seventh columns: absolute value of 
the relative error (average from 500 simulations). Eighth and ninth columns: position errors in numbers of pixels (average from 500 
simulations). 
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8 
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0.97 


1.87 


0.97 


7.97 


2.74 


(0.00,0.20;0.20) 


9 
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0.45 
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13 
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4.40 


1.13 


(0.20,0.20;0.28) 


22 
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0.13 


0.48 


0.16 


2.32 
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Table 2. NON-BLIND DETECTION. First column: pairs of values of Fq and Fu in Jy and the corresponding value of F for the 
high- noise zone (top) and the low- noise zone (bottom). Second and third columns: detection power for the NPF and FF (percentage). 
Fourth and fifth columns: flux relative errors (average from 500 simulations) for the two filters. Sixth and seventh columns: absolute 
value of the relative error (average from 500 simulations). Eighth and ninth columns: position errors in numbers of pixels (average from 
500 simulations). 
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